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Abstract. We consider a two-band superconductor with relative phase Jt between the two 
order parameters as a model for the superconducting state in ferropnictides. Within this model 
we calculate the microwave response and the NMR relaxation rate. The influence of intra- 
and interband impurity scattering beyond the Born and unitary limits is taken into account. 
We show that, depending on the scattering rate, various types of power law temperature 
dependencies of the magnetic field penetration depth and the NMR relaxation rate at low 
temperatures may take place. 
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1. Introduction 



The recent discovery of Fe-based superconducting compounds 0J has stimulated the research 
of unconventional superconductors. One of the most important and still unsettled issues is 
the symmetry of the superconducting gap function. So far, different experiments produce 
conflicting results. As regards measurements of the penetration depth and the NMR relaxation 
rate, a power law behavior at low temperatures is now clearly established, which is a 
signature of unconventional order parameter symmetry. One possible scenario of a pairing 
symmetry state is a superconductor consisting of two relatively small semimetallic Fermi 
surfaces, separated by a finite wave vector Q with the relative phase n between the two order 
parameters. This is the so-called s± model, first proposed in Ref. 0. In our previous work 
we have shown that s± model with strong impurity scattering can explain the power law 
behavior of the NMR relaxation rate. Therefore it is important to extend this formalism to 
address microwave properties of a two-band s± superconductor, in particular the magnetic 
field penetration depth and real part of complex conductivity, since experimental data are now 
available for single crystals of Fe-based superconductors. 

In this paper we calculate the microwave response and the NMR relaxation rate for a 
model s± superconductor in which impurity scattering is treated beyond the Born limit and 
discuss the relevance to the experimental data for Fe-based superconducting compounds. 



2. General expressions 

We describe a multiband superconductor in the framework of the Eliashberg approach 
equations for the renormalization function Z;(to) and complex order parameter <|>z(od). As 
shown in the first reference of [1211 . the BCS approach can give highly inaccurate results 
in the case of interband superconductivity due to the BCS neglect of mass renormalization. 
In addition there is evidence for strong-coupling in the pnictides, with many experimentally 
determined A/T c ratios substantially exceeding the BCS value of 1.76, and so we therefore 
employ the Eliashberg equations. 

On the real frequency axis they have the following form, assuming an uniform ( band- 
independent) impurity scattering (see e.g., Ref-Ol SI 13) 



4>/(o>) = I / dzKtjiz^gfiz) +i^ («t(o>) -g£(o^ 

J OO 

oo 

(Z ; (a))-l)a)=£| Jz4(z,a))Regf(z) + i^(gf(a)) + gf(Q))), (1) 

J — oo 
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where £> = 1 — a + o (gf (to) + gf (go)) + (gf (go) - gf (go)) . For our model gf (oo) = 

ni((o)Zi((o)Gi/Di((o),gf (go) =n;(Go) (|)j(Go)/Z);(Go), where D/(go) = J [Z / -(co)co] 2 - <j^(co) and 
ni (go) is a partial density of states, y = 2co/kN(0) is the normal-state scattering rate, N(0) is 
the total density of states (i.e. summed over both bands) at the Fermi level, c is the impurity 
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concentration, and a 



{■KN(0)v}- 



is the impurity strength (a — > corresponds to the Born 



l+[JiiV(0)v]' 

limit, while a = 1 to the unitary one). The kernels K^' Z (z, (fl) describe the electron-boson 
interaction and have forms 



da 



tanh + coth 
z + £1 - co - id 



where the spin-fluctuation coupling function is Bij(£l) = Bij(£l) == XijK(DQ. s f/(£l^f + Qi 2 ) 
for the equation for (]), and 2?,-y(£2) | for the equation for Z. Here X ;/ is the coupling constant 
pairing band i with band j and CIsf is the spin fluctuation frequency. Note that all retarded 
interactions enter the equations for the renormalization factor Z with a positive sign. 

We note that the implementation of the band-independent impurity scattering is contained 
in the second term on the right-hand side of Eq. 1, where the y is applied to both bands (albeit 
with a relative minus sign in the first equation due to the order parameter sign change between 
bands). We have chosen such a band-independent scattering for several reasons, including 
consistency with the previously published work and to avoid a proliferation of parameter 
choices. However, recent work of Senga and Kontani [61 suggests that this assumption is 
justified on an experimental basis. Their Fig. 4 shows that only yMer/jintra between 0.9 
and 1 is consistent with the several sets of nuclear spin relaxation rate Tj data showing 
T 2 - 5 — r 3 behavior over a very large temperature range. The theoretical rationale for such a 
comparatively large interband scattering rate remains unclear, but can be plausibly related to 
the inherent disorder in these systems, with the dopant atoms themselves acting as scattering 
centers. 

The microwave conductivity in the London (local, q = 0) limit is given by 



o'(co) = G) 2 pl / n i (co)/47c/oo, 



(2) 



where n,((o) is an analytical continuation to the real frequency axis of the polarization 
operator (see, e.g. Refs. ll7l,||l,ll9l,|fl0llfTTll) 
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and 

oo± = to' ± to/2, and the index R(A) corresponds to the retarded (advanced) brunch of the 
complex function F R ^ = ReF ± HmF ( the band index i is omitted), and (0 = Z,(co)(0. 



Here (o"f = y 87te 2 (Af,-(0)v"v^) is the plasma frequency in different directions. For the dirty 
case the low frequency limits of expressions [2] and [3] can be reduced to the strong coupling 
generalization of the famous Mattis-Bardeen expressions lfT2l 

oo 

c,((D-0) =a dc | *»(-^) {[Regf(»)] 2 + [RegfN 
o 

oo 

a dc / *»("^) {[Re^f(a))] 2 + [Regf(a))] 2 }, (4) 

where o dc = Ni(0)v 2 7 e 2 %i is a contribution to the static conductivity from i— th band. Note 
that in the London limit there are no cross-terms connected two bands. 

An important characteristic of the superconducting state is the penetration depth of the 
magnetic field X La p in the local (London) limit, which is related to the imaginary part of the 
optical conductivity by 

1 Alap = lim o 471(0 Imo ap (a), q = 0)/c 2 = (n^ReU^® = 0)/c 2 , (5) 

where a, |3 denote again Cartesian coordinates and c is the velocity of light. If we neglect 
strong-coupling effects (or, more generally, Fermi-liquid effects) then for a clean uniform 
superconductor at T = we have the relation \ L a p = c/G) pl T. Impurities and interaction effects 
drastically enhance the penetration depth, and it is suitable to introduce a so called 'superfluid 
plasma frequency' G)^ a p by the relation C0^ a p = c/X L a p. It has been often mentioned that 
this function corresponds to the charge density of the superfluid condensate, but we would 
like to point out that this is only the case for noninteracting clean systems at T = 0. 

In the two-band model we have the standard expression (neglecting vertex corrections) 

1 Ai,apOO = (< a pC0A) 2 = (6) 

v 2 

J pl,i 




<0 2 (n)+A 2 W]3/2' 



where a>(«) and A(n) are the solutions of Eq. \T\ continued to the imaginary (Matsubara) 
frequencies (A,-(«) = (j),(z(0 ;! ), G),(r) = a> n Zj(i(D n )). The calculations along these formulas can 

s f 

be thus presented in form of the effective superfluid plasma frequency, Q)^ . 

For the NMR relaxation rate, following |fT3ll , we can write down the following general 
expressions. 

Vr 1 r = 4ii mn £[F(rf^S^, (7) 

271(0^0^ (0 
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where %±(q, Cfl) is an analytical continuation to the real axis of the Fourier transform of the 
correlator 



X±(r,x) = -((r T (5 + (r,-rt)5_(0,0))), 



imp ' 



averaged over the impurity ensemble. Here S±(r,—ix) = exp(Hx)S±(r) exp(— Hx) where 
H is the electron Hamiltonian, x denotes imaginary time, and S + (r) = \|/|(r)\|/j_(r) and 
S-(r) = \|/j(r)\|/-[-(r). As a result we have 



1/TiT 
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A\k!(G>) ^;,k 2 (oo)_ 

Here ^ (co) = [coZj ; ki (co)] 2 — ^ 2 k[ — (|) 2 ki (co), ^},kiis the bare energy. For the Fermi-contact 
interaction 



1 



did 



dm 
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[Regf (to) +Regf (oo)] 2 + feegf (oo) +Regf (00 



■(9) 



This expression contains the cross-term in contrast to the microwave conductivity. In this 
paper, in the Tj calculation only these cross terms are used to emphasize the interband 
character of the superconductivity, as it is these cross terms that are most enhanced by the 
nearly antiferromagnetic state within a more detailed RPA approximation. For a single band 
system the full expression is proportional to Eq|4]when of — > oo (Ref. |fT4l ). but in multiband 
systems \/T\T and Gi (co — > 0) can behave differently. 



3. Results and discussion 



It is well known that pair-breaking impurity scattering can induce substantial sub-gap density 
of states, which can produce power-law low temperature behavior in a whole host of 
thermodynamic quantities, such as specific heat, London penetration depth, nuclear spin 
relaxation rate, and even optical conductivity. Such behavior has been well-studied in the 
two canonical limits of weak (Born) scattering and strong (unitary) scattering lfT51 , but 
the intermediate regime has received almost no attention. In addition, with the advent of 
the multiband superconductivity in MgB-2 and the apparent multiband, primarily interband 
superconductivity in the pnictides, comes a need for further study of the intermediate regime 
in an interband case. Recent studies lfT6l [T9l [T71 have addressed the effects of impurities in 
the pnictides, but only in the Born or unitary limits. Here we study the important and likely 
more realistic intermediate regime, with o, effectively the scattering strength is varied in the 



Microwave response of superconducting pnictides: extended s± scenario 



6 



O 
Q 




o i I i I i I i 

5 10 15 20 



co (meV) 



Figure 1. (color online) The quasiparticle density of states for the three indicated cases. The 
near-Born case o = 0.1 retains a small gap, while the intermediate case shows a monotonic 
DOS and the near-unitary is gapless. 



range from a = corresponding to the Born limit to o = 1 corresponding to the unitary limit. 
As stated earlier, for all calculations the impurity scattering rate y,„ /ra = y ; „ /er = 0.8Ao- 

We will now illustrate the above discussion using specific numerical models. First, we 
present numerical solutions of the Eliashberg equations using the spin-fluctuation model for 
the spectral function of the intermediate boson: 5 i; (oo) = A,;/7t(oQj//(Q^ + G) 2 ), with the 
parameters £l s f = 25meV, X\i = X22 = 0-5, and A42 = ^21 = _ 2. The rather large coupling 
constants are an attempt to model the rather large experimentally observed ratio A/T c . This 
set gives a reasonable value for T c ~ 26. 7K. A similar model was used in Ref. Il20l to describe 
optical properties of ferropnictides. This model was also used in j3l and for consistency is 
used here. As stated earlier, we further assume that each surface features the same gap []2T|. 
and that the intraband impurity scattering rate and interband scattering rate are both equal to 
O.8A0, where Ao is the low-temperature limiting value of the superconducting gap A. As in 0, 
we have chosen a relatively large impurity scattering, which is to be expected considering the 
early state of pnictide sample preparation and the limited availability of large single crystals. 

We begin with the density of states, shown below in Figure 1. Several effects are 
apparent. Firstly, for all three a values the substantial peak usually present at co = Ao (about 
6 meV here) is substantially truncated, with much spectral weight transferred below the gap. 
However, the detailed sub-gap behavior depends radically upon the scattering strength a. The 
near-Born case o = 0. 1 still retains a small minigap of approximately 1 .5 meV, which will lead 
to exponentially activated behavior below about 4 Kelvin. Although some data has shown 
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Figure 2. (color online) The inverse squared penetration depth. The near-Born limit 
approaches the BCS "two-fluid" calculation (°c 1 — T 4 ) at low temperatures, mimicking 
exponential behavior, while the other two cases show power-law behavior, as in Fig. 3. 
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Figure 3. (color online) The low temperature behavior of the penetration depth, showing the 
evolution with decreasing o from T 2 gapless behavior towards the exponential-mimicking T 4 
character. 
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evidence for such exponentially activated behavior, there is also significant data showing 
power-law behavior. The intermediate case o = 0.4 shows a monotonically increasing density 
of states and essentially no minigap, leading to power-law behavior, as proposed in J3]|. 
Finally, the near-unitary case a = 0.8 also shows a monotonically increasing density of states, 
but is nearly constant at low energy. We will see that such behavior leads to a quadratic 
temperature dependence of the penetration depth, even without the assumption of the strict 
unitary limit. Gross et. at some time ago noted [fT8l in a different context that T 2 behavior 
does not require the unitary limit. We note parenthetically that the behavior depicted depends 
rather strongly upon the large value of impurity scattering assumed; the first two cases will 
yield more exponentially activated behavior if the scattering rate is much less strong, while 
the near-unitary case can potentially [4] lead to a non-monotonic density of states. 

In Figure 2 is shown the inverted squared London penetration depth 1/X 2 (T), the 
so-called superfluid density for several cases as indicated in the figure. In all cases the 
temperature dependence of 1 fk 2 (T) is different from the standard two-fluid (Gorter-Casimir) 
model X~ 2 (r) = X~ 2 (0) 1 - (T/T c ) 4 , that is similar to the BCS result. Due to the sign 
change between gaps, the interband component of the scattering matrix is strongly pair- 
breaking, analogously to magnetic scattering in s-wave superconductors. As a result, the 
superfluid density shows near-exponential character at low temperature in the near-Born case 
(a = 0.1), while the other two cases (a = 0.4 and 0.8) exhibit power-law behavior at low T, 
with the actual power varying between 2 and 3. 

A more detailed view of the low-temperature X(T) power law behavior is presented in 
Figure 3, which shows Ak(T) /A,r=o for the same three cases. We see that the near-Born limit 
case (a = 0.1) approaches a T 4 behavior, reminiscent of a two fluid model, while the near- 
unitary case shows a fairly robust T 2 behavior and the intermediate case falls between these 
two limits, as one would naively expect. Experimental data available so far 11231 l24l I2~5l are 
consistent with either T 2 , or T 4 or exponential (gapped) behavior. Within our model, both 
results can be explained by proper choice of the impurity scattering rate. It is interesting to 
note that the T 2 dependence we obtain corresponds to strongly gapless regime. Similar results 
were obtained recently in Ref. lfT9l but in the Born limit only. 

Figure 4 shows the calculated real part of the microwave conductivity for the three cases 
above. The microwave conductivity( Fig. 4) Oi (T) does not show the coherence peak near T c . 
The suppression is connected with strong-coupling effects (see, [1221). Below T c the behavior 
of the Gi (T) is determined by the filling of the impurity induced states below A. Qualitatively 
it is similar to the temperature dependence of the NMR relaxation rate ( see Fig. 5), but in 
the latter case the Hebel-Slichter peak is additionally reduced for s± model by the different 
kind of the coherence factor. Almost all of the non-canonical BCS behavior derives from the 
interband component of the scattering matrix, which results in near constant behavior at low 
T for the near-unitary case, as might be expected from the form of Equation 51 in which a 
squared density of states enters. The intermediate case shows power law behavior as well, 
with the precise exponent not extracted. 

Finally we turn in Figure 5 to the nuclear spin relaxation rate T^ 1 for the same three o 
scenarios. Note also that following convention we have plotted (T\T)~ l rather than T^ 1 , and 
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Figure 4. (color online) The real part of the microwave conductivity. Note the substantial 
increase with scattering strength at low temperature. 




Figure 5. (color online) The temperature dependence of the relaxation rate \/T{T, exhibiting 
near-Korringa behavior for a = 0.8 and power-law behavior for the other curves. 
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all power-law references here mean (TT\)~ l . Tj has been a source of substantial controversy 
in the pnictides due to the existence of several data-sets ll26l l27l l28l |29l showing near-T 2 
behavior throughout nearly the entire temperature range, although there now exist data ll30l 
deviating from this behavior. Several things are apparent from the plot: first of all, the 
near-Born limit case shows power law behavior (\/T\T ~ T 3 ) throughout nearly the entire 
temperature range below T c , although it will ultimately revert to exponentially activated 
behavior at the lowest temperatures. Substantial impurity scattering in the Born limit can 
thus mimic much of the behavior commonly ascribed to nodes, as was noted in lfT71 [T9l . 
The intermediate case shows an approximate T 15 behavior, as was described in O, which 
is largely driven by the monotonic density of states presented in Figure 1, where the same 
parameters are chosen. Korringa behavior results in the near-unitary limit, as is again a direct 
consequence of the corresponding behavior of the density of states in Figure 1, but does not 
result in either of the first two cases unless the scattering rate y is increased significantly 
beyond 0.8Ao- 

It should now be clear that impurity scattering in various strengths (i.e, a), if sufficient 
impurity concentrations are present, can produce a wide variety of power-law behaviors in 
many thermodynamic quantities, even in the near-Born limit. In the s± state, interband 
impurities are clearly much more effective in creating such behavior. This has implications for 
the ongoing lively debate about pairing symmetry, with significant numbers of proposals for 
nodal superconductivity in the pnictides and some experimental evidence for such behavior. 

In conclusion, we have calculated the microwave response and the NMR relaxation rate 
for a superconductor in s± symmetry state by solving Eliashberg equations with a model 
spectrum and taking into account impurity scattering beyond the Born limit. We show that 
the T 2 temperature behavior of the penetration depth and the NMR relaxation rate at low 
temperatures can be reproduced in this model. We have also demonstrated the dramatic effect 
of the impurity scattering on the real part of the microwave conductivity, which in particular 
results in near constant behavior at low T for the near-unitary case. 
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